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■ Abstract 

We construct the gluon wave functions, fragmentation functions and scattering amplitudes within the light-front per- 
turbation theory. Recursion relations on the light-front are constructed for the wave functions and fragmentation 
functions, which in the latter case are the light-front analogs of the Berends-Giele recursion relations. Using general 
relations between wave functions and scattering amplitudes it is demonstrated how to obtain the maximally-helicity 
(~| ' violating amplitudes, and explicit verification of the results is based on simple examples. 
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D . 1. Introduction 



Over the past years there has been an enormous progress in the computation of the multi-particle helicity ampli- 
tudes in QCD. The tree-level gluon amplitudes turned out to have relatively simple form for special configurations 
^ \ of helicity of external on-shell gluons yj despite being the sum of a large number of Feynman diagrams. The re- 
lO ' cursion relations have been found 101 which allow the computation of amplitudes with arbitrary number of external 
\ legs. Presently there exists a number of automated tools for the evaluation of the multi-parton amplitudes for different 
O . species of particles |3i, |4, II H il ■ 

A natural question thus arises as to whether one can gain more physical insight into the structure of these results. 
Quantization on the light-front Jgt] offers particular advantages which are related to the presence of only three dynami- 
cal symmetry operators. In particular, in the infinite momentum frame there exists a remarkable isomorphism between 
a subgroup of the Poincare group and the Galilean symmetry group of non-relativistic quantum mechanics in two di- 
mensions ifloilTlll . The vacuum on the light-front is essentially structureless [IT2| (up to zero modes), which allows to 
define unambigously the partonic content of hadrons and of hadronic wave functions. The simplicity of the vacuum 
\ has been used to argue about the presence of in-hadron quark condensatesJ131]. The light- front framework has been 
H ■ used to investigate the hadron dynamics from AdS/CFT correspondence lll4ll . Methods of the light-front perturbation 
. . .' theory (LFPT) were used to compute the soft gluon component of the heavy onium wave function and to obtain a cor- 



respondence with the hard Pomeron in QCD 111511 . Similarly, the LFPT has been utilized in the small x limit to analyze 



the quark distribution function of a heavy meson utilizing ladder operators |16| |. In work 111711 light-front perturbation 
techniques yj, 19, 20, 21 J have been applied to construct the Hght-front wave function of the gluon. Expressions for 
the recursion relations which are obeyed by these wave functions have been found. They have been solved for some 
special cases of the helicity configurations. To be precise, an exact form of the tree-level multi gluon components of 
the real gluon wave function has been found for special helicity configuration in which both the initial gluon and the 
final gluons carry the same helicities. It turns out that the variables which occur naturally in the light-front formulation 
are closely related to the spinor products which are used to express the results for the maximally-helicity-violating 
(MHV) amplitudes fll E2|] . The scattering amplitude for 2 — > n process with massless gluons has been also considered 
in the LFPT, under the special assumption that the evolved projectile gluon is separated in rapidity from the target 



gluon lllTll . The result in this approximation has been shown to coincide with the expression of the Parke -Taylor 
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amplitude. As an additional result the closed form for the fragmentation function of an off-shell gluon into arbitrary 
number of final state gluons have been computed. A recursion relation for the fragmentation function has also been 
found. The recursion involves factorization of the fragmentation trees which is not apparent in the LCPT due to the 
entanglement of the particle momenta in the energy denominators for the intermediate states. It has been shown that 
the factorization of the fragmentation trees is achieved after the summation of all the light-front time-orderings of all 
the splittings in the cascade. 

The main goal of this paper is to rederive, using the LFPT and techniques presented in lll7ll for fragmentation 
functions and wave functions, the exact expressions for the helicity multi-gluon amplitudes in QCD. In particular it is 
important to relax the simplifying assumption used in Ref. IitIi about the rapidity difference between the projectile and 
the target gluon. We first argue that the factorization property of the fragmentation functions is the light-front analog 
of the Berends-Giele recursion relations for the gluonic currents. A generalized factorization for the fragmentation 
function is constructed which includes more complicated helicity configurations. Next, the explicit relation between 
the light front wave functions and scattering amplitudes is shown and verified by the computation of the lowest order 
examples. Additionally it is shown how the recursion relations for the wave functions in the case of the same helicities 
for the outgoing gluons are related to the vanishing of the scattering amplitudes with -,+,...,+ configurations. The 
vanishing of these amplitudes is related to the kinematical constraint of the conservation of the light-front energy 
between the initial and final state. The vanishing of the amplitudes with the same helicities for all the gluons is related 
to the explicit conservation of the total angular momentum on the light-front. Finally, the recursion relation for the 
wave functions in the case when one of the outgoing gluons has different helicity is constructed which leads to the 
recursion relation for the first nontrivial set of MHV amplitudes. Through the explicit computation the Parke-Taylor 
amplitude is reproduced. 

The structure of the paper is the following; in the next section we briefly summarize the results and techniques 
derived in the paper jlln . It is then argued that the factorization for the fragmentation functions is the light-front 
analog of the Berends-Giele recursion formula. In Sec. 3 the general relations between the light-front wave function 
and the scattering amplitudes are discussed and examples which confirm this identification are shown. The vanishing 
of the +,..., + and -,+,..., + amplitudes on the light-front is demonstrated. The recursion relation for wave function 
is presented and it is shown that it can be used to rederive the Parke-Taylor amplitude. Finally, in the last section we 
state conclusions and the outlook for future work. 



2. Gluon wave function and fragmentation function with exact kinematics 

In this section we shall first set up the notation which will be used through the rest of the paper Next, we will 
summarize the most important results which were derived in paper UtIi . In particular, we shall show the exact form 
of the «-gluon components of the gluon wave function, ^F,,, « = 1,2, .. . and fragmentation function T„ using the 
LCPT for a special choice of helicities. In the following by the gluon wave function we mean a sub-amplitude on 
the light-front where the incoming gluon is taken on-shell and the outgoing gluons are all off-shell. Similarly, by 
the fragmentation function we mean a sub-amplitude with the initial gluon off'-shell and the outgoing gluons on-shell. 
Both wave functions and fragmentation functions can be used as building blocks to construct the scattering amplitudes 
and they obey separate recursion relations. 

For now, we shall focus first on the gluon wave functions and as in Ref. ITtIi we assume that the incoming on-shell 
gluon has a positive helicity and all the final off-shell gluons in the wave function have all positive helicities. The 
resulting resummed gluon wave function will have a very simple form. 

We consider an incoming gluon with the four momentum P, and the color index a, which evolves into a state 
containing n gluons with momenta {ki, . . . k„) and color indices (ai, 02, . . . , a„) correspondingly. This is illustrated in 
Fig- El The momentum for the initial virtual gluon with virtu ality -Q^ and a transverse momentum q can be written 
in the light-front variables, 

with P"^ - P'-' + P^ and zo - ^- The momenta of the last n - gluons are labeled ki, . . . ,k„, as in Fig [T] Each 
of this momenta can be represented as ^ = {ziP^ ,kT ,k.), with Zi being the fraction of the initial P^ momentum 
which is carried by the gluon labeled by ; and k. being the transverse component of the gluon momentum. In the 
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Figure 1 : The multi-gluon wave function of the incoming gluon with momentum P. The vertical dashed lines sym- 
bolize different intermediate states where we need to evaluate the energy denominators. It is understood that the wave 
function scatters finally on some target. 



LFPT dill, Sin one has to evaluate the energy denominators for each of the intermediate states for the process 
depicted in Fig. [T] The energy denominator for the n gluons in the intermediate state is defined to be a difference 
between the light-front energies of the initial and the intermediate state in question. For the wave function shown in 
Fig. [T] we assume the last intermediate state is with n gluons. The corresponding energy denominator for this state 
reads 



D„ = P 



" 1 
/_i I p+ 



— - — + — + — 4-... + — 



Zl Z2 



where the light-front energy of the gluon / is equal to 



D„ 

p+ " 
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and we introduced the auxiliary notation for the (rescaled) denominator D„. 

The wave function shown in example in Fig. [1] would thus be given schematically by the expression 



where Vj are the vertices and zj and Dj are the corresponding fractional momenta and denominators for all the 
intermediate states. 



The results derived in 111711 and in the following sections are for the color ordered multi-gluon amplitudes. There- 
fore we focus only on the kinematical parts of the subamplitudes. 

In Ref. d only a special component of the gluon wave function was considered which lead to many simplifica- 
tions and allowed for the recursion formulae to be solved relatively easily. Namely, it was assumed that the gluons 
have all the same (e.g. positive) helicites. The convention used in this work is such that the leftmost gluons have 
incoming momenta and rightmost gluons have outgoing momenta and the light-front time on figures flows from left to 
right. Following lll7ll we choose to work in the light-cone gauge, 77 ■ A = 0, with vector rj^^ - (0, 2, 0) in the light-front 
coordinates. The polarization four-vectors of the gluon with four-momentum k have the following form in this gauge 
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where - (0, 0,6'"^^), and the transverse vector is defined by - ±0- Due to the intrinsic properties of 

the three and four gluon vertices when projected onto heHcity states, the resuhing gluon wave function with the same 
(positive) gluon heHcities in the final step can be formed with the gluons with positive helicities at all intermediate 
steps. The 4-gluon vertex vanishes for this helicity configuration and therefore the splittings are given only by the 
3-gluon vertex with a (H — > ++) helicity projection, which corresponds to the situation depicted in Fig. [T] For the 
case of interest the projection of the vertex reads 

V%Tikuk2,k3) = 2g6{zi+Z2+Z3)6^\k^+k,+k,)r"'-"'zi ■ - ^) . (3) 

~ ~- " \Z3 Zl) 

Here g is the strong coupling constant and f"'"^-"^ is the structure constant for the SU{Nc) color group. In the 
above notation we assumed all gluons to be outgoing, therefore the gluon 1 has a negative helicity which translates 
onto positive one when inverting the momentum flow. For the case of interest, namely that of (h — » ++) transition, the 
amplitude is described by Eq. Q, with z\ being the fraction of the + component of the momentum of the incoming 
gluon. Then the vertex is proportional to the variable 

(4) 

\Z3 Zlj 

It is well known lflol[l8l[l9ll that on on the light-front the Poincare group can be decomposed onto a subgroup which 
contains the Galilean-like nonrelativistic dynamics in 2-dimensions. The '+' components of the momenta can be 
interpreted as the 'masses'. In this case the variable (2) can be interpreted as a relative transverse light-front velocity 
of the two gluons. Interestingly enough, the same variable is present when we consider the change of the energy 
denominator due to the splitting. The importance of this variable is evident when considering the explicit forms of the 
energy denominators. For example for the situation depicted in Fig. [T]the difference in the energy denominators for 
the last and last-but-one intermediate state can be shown to be equal 

D„ - D„_i = = ^4 . (5) 

Here we introduced another combined variable 

- ^ , (6) 

Zi + Zj 

which de pen ds only on the longitudinal degrees of freedom. It can be interpreted as the reduced mass for a two body 
problem (isl]. Thus the difference in the energy denominators before and after the splitting is proportional to the 
(twice) kinetic energy associated with the relative motion of particles 3 and 4. 

In the previous work we have adopted the convention where the normalization factors for the intermediate lines 
were absorbed into the vertices, i.e. 

VTT7(kuk2,ki) = , V"'ff(kuk2,ko,) = 2gf'"'-"' . (7) 

Using the definitions of the vertex and the energy denominators presented above one can form the recurrence 
formula for the gluon wave function with n components 

lie e*^'v.. , 

%(^i, . . . , ki, ki^u- ■ ■ , k„) = -p^ ^ - , %-i(^i,. . . , kiM,. . . , k„) . (8) 

Here, '^n-\{k\,k2, . . . , . . . , A:„) (with ka+x = kj+kj+i) is the (n-l)-gluon wave function in momentum space before 
the splitting of gluon with momentum , and ^„{ki ,k2,..., ki, kj+x, . . . , k„) is the wave function after splitting of 
this gluon. To obtain the full wave function one needs to sum over the different possibilities of the splittings which 
gives the following result 



where D„ is the denominator for the last of all the intermediate states with n gluons. We note that all dependence 
on momenta of daughter gluons / and / - 1 is now encoded in two variables: , and v^. j .. The recurrence formula 
will be utilized for the derivation of the scattering amplitudes. We note that this recurrence formula is written for the 
special choice of the helicities. It can be readily generalized for the other case of hehcities. We will come back to this 
point later when we use it to derive the Parke-Taylor amplitude. 

The variables v^.^ which represent the 'velocities' of gluons on the light-front are actually related to the variables 
used in the framework of helicity amplitudes, see ll22ll for a nice review. For a given pair of on-shell momenta A:,- and 



,k k.\ Ik k 

{ij)- yf^ji-^ ■{---], V^l^"'- ---h (10) 



kj we have the result 



where the latter are defined by 

= + (11) 
The chiral projections of the spinors for massless particles are defined as 

|/+) = ^^{kd = i(l±75)<A(^,) , (±/| = ^pJkd, (12) 

for a given momentum kj. The spinor products are complex square roots of the total energy mass squared for the pair 
of gluons (/, j) 

{ij)[ij]-{ki + kj)\ (13) 
and they also satisfy {ij) - {ij}*. Using the above definitions (110113b we have that 



IK.: ^1,'^ 



k. k. 



,Zi Zjl 

which is real and positive for the on-shell gluon momenta. Finally, combining (|4]i and ( fTOt we obtain 

< ij > = ■ V,., , [ U ] = e*^' ■ V,,. , (14) 



and the dependence on the transverse momenta in the light-front wave function can be expressed by spinor products 
<;;) and [/;■]. 

2.1. The multi-gluon wave function in the case of the on-shell incoming gluon 

The multi-gluon wave function *F„ has a closed compact form in the case of an on-shell incoming gluon, ^ 
with helicity + and when all the final gluons also have + helicities. The simplification stems from the fact that the 
energy denominators do not contain the virtuality Q^. 

The recurrence formula, i.e. Eq. (|9), for the case of the on-shell incoming gluon and the same helicities for the 
outgoing gluons has the explicit form 

A,+i*l'«+i(l,2,...,« + l) = 
= 2/g^>P„(12,3,...,«+l) + 2ig^'V„(\,23,...,n + \) + ... + 2/^^|^»I'„(l, 2, . . . , nn + 1) , (15) 

with D„+i - l^i/z\ + k^lzi + ■ ■ ■ + /Zn+i - Izi). We have introduced the notation *P„(1, ...,/- !/,...,«+ 1) 
where / - 1 / means that it is the gluon with the transverse momentum ; = + ki. Again, we recall that we are 
considering color ordering in the amplitudes, therefore we will suppress color degrees of freedom. The wave function 
for the incoming state has the normalization ^i(l) = 1. After the first splitting one gets 

^2(1,2) = 2ig-^^ = -igJ^-!—, (16) 

V^12 ^12Vi2 Vfl2 ?12V12 
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(a) 



(b) 



Figure 2: Graphs contributing to the gluon wave function with 3 gluons. 



where we have taken that Qp' - 0. According to ( fTsT i, the next splitting leads from *F2(1,2) to *F3(1,2,3) with the 
graphs depicted in Fig.|2]and the result is 



£)3*P3(1,2,3) = lig 

1 V* 



1 



Vfl2^(12)3 ^( 1 2)3 V(l 2)3 V^23^H23) ^ 1 (23) V i (23) 



*P2(12,3) + 



"23 



^2(1,23) 



VI23 

V^O V'*2^1(23)Vl(23) + V23^(12)3V(12)3 



(17) 



where we used ^i2^(i2)3 = fe^i(23) = 
obtain. 



^3(1,2,3) - -g 



V^TE^I ^(12)3^1(23)^(12)3 Vi(23) 

- £i±2±3 can further simplify the numerator of this expression to 
Vz^y 1 1 



(18) 



y]Z\Z2Zi ^(12)3^1(23) l'(12)3 l'l(23) 

Note that, the energy denominator D3 disappeared from the equation as it has canceled with the numerator when 
finding the common denominator for expression (ITTt . This is actually quite an important property and leads to 
the very simple expression for the resummed gluon wave function. It is also important in recovering some crucial 
properties of the helicity amplitudes as we shall show later. The resummed wave function with n-gluon components 
can be showed to be given by the following compact expression IPtIi 



>F„(1,2,...,«) 



(-l)"-'(/g)"-'A^"> 



1 



^^!z\zJl~^^ ^(12...n-l)« ^(12...n-2)(n-l n) 
1 



^1(2...«) 



(19) 

V(12...«-l)« V(12...«-2)(«-l «) ■ ■ • V\(l...n) 

Throughout this section we have been using the following complex representation of the transverse vectors: v,^ - 

^/] + + - . • + kj^ kj^ + . . . + kj 



e* ^ ■ V. ., V*. = e^^' - v ., and a useful notation. 



Zl +Zi,+ ...+ Z: 



Zu +Zi. + . 



+ Z 



7, 



{Zi,+Zi, + ... + Zi;}{Zj,+Zp_ + ...+Zj,) 



Zi, + Zi, + . . . + Z,„ + Zj, + Zj^ 

together with ki = e*^* ■ k_.. We also introduced notation for the partial sums Z{\...i) = Z\ + Zi + ■ ■ ■ + Zi and 
= ki + k2 + ■ ■ ■ + kj. In formula ( fT9] l, A'"' is defined to contain the global momentum conservation 5-functions, 
^(n) = ^(2) ^ k. - ^ ^ ( 1 ~ 2"=i Zi ). Recall, that on the light-front, the longitudinal momenta and the transverse 
momenta are conserved but not the light-front energies k^ (except for final and initial states). 

The formula (fT9l l can be proven by induction (for details see iITtIi ') and with the help of the following identity 



(20) 



(21) 



V(l.../)(i+l...n)^(l...i)(i+l 



7=1 



2(l...n) 



7=1 



(1-0 ■ 



^(1-0 

Z(l...n) 



kn...n) - ki 



(1.../) 



^(1-0 
( 

20 



(22) 



where the momentum conservation conditions q - k^i,,,,,) and zq - Z(i...„) have been utiHzed. The above identity can be 
used to show that 

-\d, ■ (23) 



2^ l'/,+ l«I.../)(/+l...n)^(l...i)(i+l...n) - -T 

,=1 ^ v,=i 



The above expression is what appears in the numerator after combining the different terms in the recursion formula 
to a common denominator It therefore cancels with an overall denominator for the wave function leading to the simple 
form of the wave function (fT9l l. The example of that property was used to simplify Eq. ( fTTT i. 

2.2. Fragmentation of an off-shell gluon 

Using the techniques developed for the wave functions one can similarly resum the graphs in a kinematical situ- 
ation when the initial state is an off-shell gluon, and then it fragments to give the final «-gluons which are on-shell. 
There are interesting similarities between the two cases, i.e. wave functions and fragmentation functions, which can 
be expected as the topology of the graphs is exactly the same. There are some important differences however, which 
are related to the properties of factorization. Namely, the fragmentation functions can be shown to factorize whereas 
the gluon wave functions do not exhibit this property (at least not in momentum space). Using the notation of the 
previous section one can show that for the case of the splitting of an off-shell gluon, denoted by (12), into 2 on-shell 
gluons (denoted respectively by 1 and 2) the fragmentation amplitude can be expressed as 



2'g V\2 _ ^12 _ . £-3/2 J_ 



r2[(12) ^1,2] = -^-^ = -2ig = igf,^'^ — , (24) 



where D2 = Zi?=i - (^(12)) /^(12) for this case. The general expression for the fragmentation part of the amplitude 
for 1 to n gluons as r„[(12 . . . «) ^ 1, 2, . . . , n] on a tree-level is given by lll7ll 



\3/2 , 

r„[(12...«)^l,2,...,«] ^ (/g)"-'( . (25) 

\ZlZ2---Z„ V12V23 . . . V„_i„ 



The off-shell gluon is labeled by (12 . . .n) and the final state n on-shell gluons are labeled as 1,2, . . .n. These n 
final state gluons have transverse momenta . . . and the longitudinal momentum fractions z\,- ■ ■,z„- The initial 
gluon has transverse momentum k^^ and the longitudinal fraction Z(\...„) where we again used shortcut notation 
^1 „) - 2'j=i Kj ^^'^ Z(i2...n) = Zi/Li Zi ■ The fragmentation function is depicted in Fig. |3] 




Figure 3: Pictorial representation of the fragmentation amplitude r[(12 . . . «) — > 1, 2, . . . , «] for the single off-shell 
initial gluon. 

The above formula can be again proven by mathematical induction using analogous relations to the (l22l i and 
(|23]) . On top of that, one crucial assumption is that of the factorization property of the fragmentation functions 
which are formed from topologically disconnected trees and originate from different off-shell parents. Physically, the 
factorization property of independent tree amplitudes should be quite intuitive. One considers a fragmentation tree 
with a topology = ©i U . . . U 0,„ where ©i, . . . &„, are topologies of the fragmentation trees of m parent virtual 
gluons into n on-shell gluons, so that the first parent gluon, denoted by (12 . . . ni) fragments into gluons (1,2,... , «i), 
the second parent gluon, («i + 1, «i + 2, . . . , ni) fragments into gluons («i + 1, «i + 2, . . . , «2), and so on. One can then 
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show that the amplitude of the fragmentation tree 0, denoted by Tq may be factorized into fragmentation amplitudes 
of the parent gluons, Tq^ in the following way, Tq[{1 . . . («i + 1 . . .«2), ■ ■ ■ , («m-i • ■ • «) — > 1,2, .. .n]. 



r0[(l . . .ni),(ni + 1 . . .«2), 
X r0j(«i + 1 ...ni) 



,(«„,_!...«) -> 1,2, ...n] = r0,[(l...ni) ^ 1,2,. 

ni + l,ni + 2, . . .«2] X . . . X Te [(«„_! + 1 . . .n) ^ nm-i,nm-i + 1, 



(26) 



The factorization property is trivial in the case of standard time-ordered perturbation theory. However, it is not 
obvious in the light-front formulation, as for a given time-ordering of the splittings, the variables related to different 
trees are correlated in the energy denominators, and the factorization property does not hold for an individual diagram 
with particular ordering of the vertices on the light-front. It only holds if the complete summation over all possible 
orderings of the light-front time-orderings is performed within all connected trees, while preserving the topologies as 
has been demonstrated in lIlTIl . 

The above factorization property can be used then to write down the explicit recursion formula for the fragmenta- 
tion amplitudes. Namely, the fragmentation into n+l gluons denoted by T„+i [(1,2, . . . , n + 1) — > 1, 2, . . . , n H- 1] can 
be represented by lower fragmentation factors r/[(l . . . /) 1, . . . , / ] and r„+i_/[(/ + 1 . . .n + 1) — > i + 1, . . . , n -I- 1] 
and by summing over the splitting combinations. To be precise one has 



r„+i[(i2...« + i)^ 1,2,. 



)(/+l...n+l) 



X r,[(i ...0^1,...,/] r„+i_,[(/ + 1 . . . « + 1) ^ / + 1, 



,n + 1] 



(27) 



This expression is the final state analog of formula (|9]l for the iteration of the wave function and it is schematically 
depicted in Fig.|4] 



fc(12...n + l) ' ^' 




Figure 4: Pictorial representation of the factorization property represented in Eq. (l27t . a light-front analog of the 
Berends-Giele recursion relation. The helicities of the outgoing gluons are chosen to be the same in this particular 
case. Dashed vertical line indicates the energy denominator D„+i. 



The above defined fragmentation functions T can be related to the gluonic currents which play an important role 
in the Berends-Giele recursion relations. These recursive relations utilize the (gauge-dependent) current J'', which 
is obtained from the amplitudes by taking one of the gluons off-shell. The dual subamplitudes can be obtained by 
contraction with the polarization vector and putting the gluon back on-shell, ll22i1 

M(0, 1, 2, ...,«) = iP^e^Jf,{l,2, n)p=^p, (28) 

where we have defined P - Pn and in this formula p, denote the four- vectors for the momenta of the outgo- 
ing particles. In the case of the light-front calculation the light-front current can also be defined and related to the 
fragmentation function in the following way 



r„((l , 2, . . . , n) -» 1 , 2, . . . , n) s e^'(12 . . . «)^(1 , 2,...,n) 
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where by e'^(12 . . .«) we denote the polarization vector of the incoming (off-shell) gluon in the fragmentation func- 
tion. With such definition the factorization property for the fragmentation function ( ITTl i is a light-front analog of the 
Berends-Giele 101 recursion formula (see also [|23i l for the recurrence relations in the light-cone gauge) 

. n-l 

J^iX,!, ...,n)^~Yj K\pi-i^PM-nW,- ■ ■ , + 1, . . . , «) 

/=1 

n—l n-2 

'=./+I J=l 

The simpler form of (IZTT i (as compared to (|29}) which only includes three gluon vertex stems from the fact that it has 
been written for particular configuration of helicities. It is possible to write down a general factorization (recursion) 
relation for the fragmentation function which will include the 4-gluon vertex as well as the Coulomb term. For 
example in the case of the fragmentation into n-\ gluons with positive helicity and 1 gluon with negative helicity the 
recursion relation of the fragmentation function has the following form 

T ( " V* 

-r r/n , i\+ 1- o+ ,1+1 ^'■S {l...i)(i+l...n+l) 

r„+i[(12 . . .« + 1) — » 1 ,2 ,...,«+! ] - 



A,+i \ii 



~i V^(l-')('+l-n+l) 

X r,[(l .../)+ ^ r,2+,...,/+] r„+i_i[(iH- 1 ...n-H 1)+ i+l^,...,n+l^] 



n 



+ V v(i...,xi...«+i) J r,[(i . . . 0" -> r, 2+, . . . , r ] r„+i_,[(/ + 1 ...„ + 1)+ ^ / + 1+, ...,„ + 1+] ] 



;=1 



(30) 

In the above formula the vertex V4 contains both 4-gluon vertex and the Coulomb term 

,,-++_ 2;^^ Zi+l...j - Zl...i Zj+l...n+l 

^jZ\...n+l Zl...i Zi+l...j Zj+\...n+\ Zj j 

As mentioned above the Berends-Giele relations can be written on the level of individual diagrams, whereas for the 
derivation of the analogous recursion relations on the light-front (l27l i. ( l30t the summation over the time-ordering is 
necessary to decouple the disconnected fragmentation trees. 

Thus using wave functions and fragmentation functions one can write two distinct recursion relation in the LFPT 
There is however an important difference between the recurrence formulae in the two cases on the light-front, i.e. 
fragmentation and wave functions. The factorization property enjoyed by the fragmentation amplitudes does not hold 
in the case of the wave functions. This stems from the fact that in the latter case one still has the energy denominator 
corresponding to the n gluons, see Fig.[I] This means that the n gluons of the wave function are in the intermediate state 
and still have to interact. Therefore they do not necessarily form disconnected trees and therefore can be correlated 
and cannot be factorized. Still, the recursion formulae for the gluon wave functions will prove to be very useful, and 
can be used to obtain the helicity amplitudes as we will see in the following sections. 



3. MHV scattering amplitudes from the light-front wave function 

In the previous sections we have analyzed the (tree-level) wave function with arbitrary number of gluons in the 
light-front formulation as well as the fragmentation function of the off-shell gluon into n final state gluons. The 
main goal of this section is to show how to reproduce on the light-front the helicity amplitudes from the previously 
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derived wave functions. We recall that the Parke-Taylor amplitudes 111] (see |22| for a comprehensive review) have 
the following form 

M, = 2 tr(f'"f°^..f«")»j(pi,ei;/?2,e2;...;p„,e„), (32) 

where 01,02, .. o,„ p\,p2, ■ ■ Pn and ei, 62, . . . , e„ are the color indices, momenta and the helicities of the external 
n gluons. Matrices f are in the fundamental representation of the color S U(Nc) group. The sum in (l32l l is over the 
{n - 1)! non-cyclic permutations of the set {0, 1, ... , «). 

The kinematical parts of the amplitude denoted by m(l, 2, . . .,n) = m{p\, e\; p2, 62; • ■ • ; Pn, e«) are color inde- 
pendent and gauge invariant and enjoy a number of important properties. Among them is the requirement that the 
amplitudes where all the gluons have the same helicities, or only one is different from the others are vanishing (here 
we used the convention that all the gluons are outgoing) 

m(±, +,..., ±) = m(+, +, ±, ...,+) = . 



These properties can be demonstrated on the tree-level using the supersymmetry relations Il24ll . The non-vanishing 
amplitude for the configuration (-, -, -H, . . . , -h) at the tree-level is given by the formula 



m(r,2-,3+,...,«+) = ig"-^ — , (33) 

* (12)(23>...<n-2n-l)<n-ln)<nl> 



where the spinor products are defined by Eqs. (ITTl . dTTi 



In the previous work H17I1 a simplified derivation of the Parke-Taylor amplitudes on the light-front has been demon- 
strated. The simplification originated from the fact that the projectile gluon with the evolved wave function was al- 
lowed to scatter off" the target gluon which was separated by large rapidity. This allowed the exchange between the 
projectile and the target to be be treated in the high energy limit. In this limit the interaction between the projectile 
and the target is mediated by an instantaneous part of the gluon propagator in the light-cone gauge defined as 

inst (^+)2 ■ 

The kinematics of the exchange was simplified but the internal structure of the projectile was treated without any 
simplifications. In addition after the interaction the evolution in the final state was also included. To be precise the 
fragmentation of the gluons after the instantaneous interaction took place was included through the fragmentation 
functions T defined above. 

After taking this simplification it can be shown that the Parke -Taylor amplitude can be represented in this limit in 
the factorized form 

Mia,Q-^ b,l,2,...n) = -^=g'- t„(l,2, . . .;i). (34) 

The factor ~ g^s/\t\ comes form the gluon exchange and *P„ is the amplitude of the gluon evolution of the projectile 
gluon into «-gluon final state. The factor 1 / y/ZaZt is present because of the convention. The expression for the general 
n for the subamplitude can be found from the following formula 

n 

t„(l,2,. ..,«) = ^ ^ >l'„,((l...ni)(ni + l...n2)...(«m-i + !...«)) 

m—l ( !<«! <n2<...<n,„_i <h) 

X r[(l ...«i) ^ r[(«i + 1 ...«2) ^ni + l,...,n2] ... T[{n,„^i H- 1 ...«) ^ «,„_i + l,...,n] . (35) 

This formula is obtained by taking into account all possible attachments of the exchanged gluon in the cascade. The 
additional assumption is the factorization of fragmentation of virtual gluon ( l26l l. 

The final result is the following general form for the ^i*,, for an arbitrary number of emitted gluons 

%(l,2,...n) = ^"-1^^^=^^ 1 . (36) 

kl/Zl ^/ZlZ2 ■■■Zn Z1Z2 ■■■Z„ Vi2V23...V„_i„ 
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which can be shown to arise lll7ll when the recursion formula (1351 ) is solved using the explicit expressions for the 
fragmentation functions and the wave functions. Finally it was shown that the (approximated) expression for the 
Parke-Taylor amplitude can be re-derived in this Umit by combining Eqs. ( l34b and ( |36] ). 

However, as mentioned above this derivation is only approximate since it is based on the assumption of the large 
rapidity difference between the incoming projectile and the target gluon, which is only appropriate in the high- energy 
limit. In the following, we shall show a different method which will relate the wave functions with the amplitudes and 
allow to obtain the exact results without any kinematical approximations. 



3.1. Simple example: 2^2 amplitude on the light-front 

In order to prepare the ground for the derivation of general relations between gluon wave functions derived earlier 
and scattering amplitudes, we shall consider at first a lowest order example of the 2 — > 2 amplitude using the rules of 
the LFPT and the notation introduced in the previous sections. The goal of this simple example is to illustrate the fact 
that the resulting amplitude has a very similar form to the gluon wave function for 1 — > 3 transition. 




(d) 



(e) 



Figure 5: Basic graphs which contribute to the 2 — > 2 amplitude in LFPT. Dashed lines indicate the energy denom- 
inators for the intermediate states. The cross on the vertical line in graph (e) indicates the instantaneous - Coulomb 
gluon exchange. 

To calculate the sub-amplitude for 2 — > 2 gluon scattering M2^2 one needs to compute and add up all the graphs 
shown in Fig.|5] Let us focus on the situation where the helicities for this amplitude are m(0^, 1,2^, 3^), the momenta 
of the gluons and 1 are incoming and the momenta of the gluons 2 and 3 are outgoing. It is well known that the 
amplitude for such helicity configuration has to vanish, see 1221 and discussion in the previous section. Thus it is 
the simplest example, nevertheless it is instructive to see how that happens on the light-front, using the notation 
which was given in the previous section. Graphs d and e give contributions which vanish identically for this helicity 
configuration. However, the contributions from graphs a, b, c do not vanish separately. The sub-amplitude Ga for the 
graph in Fia. |5a]is given by 

/ V1V2 

(37) 



Z21 D3 



where the vertices and the energy denominator are given by 















(-- 







3(21) 



(38) 
(39) 
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D3 

to give 



§0 _ §(21) _ S 
Z(2l) Z3 



(40) 



n /I - 2 3(21) 21 

Ga = 4;g zo- ;i ■ (41) 

where we used notation 1 to denote the change in the sign of the corresponding transverse and longitudinal momenta 

k k -k 

i.e. ZQ = Zi - Zj and ^.jj = k. - kj and V3 (2I) - ^ - ^^T^' ^1^° have dropped an overall, dimensional factor of P'^ 
from the energy denominator to facilitate the comparison with the previous expressions. Similarly, the sub-amplitudes 
Gb and Gc for the graphs in Figs. |5b]and|5c] respectively, are given by 

^ A - 2 ~^12)3^21 

Gb = 4ig ZO- , (42) 

Zl ~(l2t Z2 

„ ..2 ^1(23)^32 

= 4!g zo-T — -r- ■ (43) 

io _|_ S_ _ ^23) 
^0 Zl Z(23) 

We have to take into account that the internal lines need to have the positive fractions of the longitudinal momenta in 
Figs. |5a]and|5b] This means we need to include the corresponding step functions which enforce the ordering of the 
momenta. Hence, M2^2 resulting from these three graphs is given by 

M2^2 = @iZO - Z3)G, + 0(Z3 - Zo)Gb + Gc . (44) 

It is trivial to show that Ga and Gb are actually the same and can be combined into one expression. Using 
^ + £- k + ^vve can write the denominator in (l42t as -— - — + — . The following identities also hold 

zq Zi Z2 Z3 Co C(|2] ;:3 ^ 

^(12)3 = ~^3(2i)' -2(12) = ~^(2\) ■ (45) 

Using them in ( |4TI ) and ( l42l i we realize that the expressions for Ga and Gb are the same which gives 

M2^2 = Go + G, , (46) 

where the theta functions are absent (only global momentum conservation conditions on external longitudinal mo- 
menta are present but they are implicit). The fact that the graphs a and b give the same contributions, can be shown 
to hold for more complicated graphs. This is discussed in detail in Appendix A and it will be used later to prove the 
relationships between wave functions and scattering amplitudes. 

Having reduced the expression for the sum of both graphs to (|46] | we proceed to show that amplitude M2^2 in ( l46] l 
has identical structure to the form of the gluon wave function *i'3 in Eq. ( fTTI i (modulo overall factors). To demonstrate 
this one can use the relation D3 = ^3(2i)li!3(2i)l~ ^'^ obtain 

G„ = -2/^2^0^^^^^ , (47) 

S3 (21)'^3 (21) 

Similarly, we can rewrite the energy denominator of graph (c) 

^ —1 ^23) 2 
' ~ ~^(32)I IZ(32)I I ' 



to get 



G, = -2ig'zo—^ . (48) 

S(32)I^(32)I 
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Putting both expressions together we obtain. 



M2^2 = -2ig'zo + (49) 

\?3 (21)^3(21) S(32)l^(32)l / 
T z / ^'21^(32)1 ^(32)1 + ^32^3 (21)^'3 (21) ^ 

-2i8 zo I . (50) 

S3 (2l)S(32)P'3 (21)^'(32)I / 



We notice that the structure of M2^2 now has exactly the same structure as the wave function for in (fTTb . when 
replacing 1 — » 1 . The only difference is the overall normalization factor due to slightly different convention used in 



OlTIl to derive ( fTTl i. Most importantly, just as in (fTTT i. the numerator in (l50t is proportional to the difference of the 
energies from the initial to the last state in the problem. Since in the case considered (i.e. on-shell amplitude) the last 
state is on-shell (final) state, the numerator, and thus Mo^i, will be zero. Therefore we see that the relatively compact 
form of the gluon wave functions derived in |l7| and the recursion relations for these wave functions are closely 
related to the properties of the amplitudes. In this particular case the simplicity of the wave functions are related to 
the fact that the amplitudes with a single helicity different from the rest are vanishing. The result derived for lowest 
order scattering can be generalized to arbitrary number of the final state particles (at the tree-level). It is important to 
note though that the wave functions derived in Sec. |2]themselves are non-zero due to the off-shellness of the last state. 

3.2. General relations between the wave functions and the amplitudes 

Given the result of the explicit example shown in the previous section we proceed to investigate the relation 
between the amplitude for a general number of final state gluons M2^„ and the wave function ^„+\ which were 



derived previously, see Ref. lll7ll and Sec. |2] That such relation should exist is quite intuitive as the two objects 
correspond to the graphs which are topologically equivalent with only the one external leg changing the momentum 
from the final to initial state. In the case of the LFPT it is however not an entirely obvious relation as the energy 
denominators, which entangle the momenta of the particles in a given graph, in principle change via such operation 
and, moreover, some graphs are vanishing. For example, for the case of the *P3 we have only two graphs depicted in 
Fig.|2]which contain the 3-gluon vertices and give non-trivial contributions, whereas for the 2 — > 2 amplitude there are 
three graphs (with 3-gluon vertices) whose contributions are non-vanishing, those shown in Fig.[5al|5b]|5c] Therefore 
by changing the momentum line from outgoing to incoming in Fig. |2] one needs to perform additional summation 
over the time-ordering of some of the vertices, so that in our example the contribution from graph (a) in Fig. |2]has to 
give two separate contributions in Fig. |5a]and|5b] We shall show how to identify and sum over the time-ordering for 
the case of 2 — > 3 and then we perform a proof for the general case and show that the amplitude from the M2^„ graphs 
can be obtained through analytical continuation from the fight-front wave functions ^n+\ which contain in general a 
smaller number of graphs. 

We will restrict ourselves to the case of only 3-gluon vertices which will limit the possible helicity states. However, 
the methods developed are much more general and we will utilize them later to find amplitudes with more complicated 
heUcity states. 

Let us first consider the case of 2 to 3 gluon scattering. We want to find the M2^3 graphs from the graphs of the 
gluon wave function which describes 1 to 4 gluon transition in the H — > + + ■■■+ configuration. This gluon wave 
function contains (4 - 1)! =6 graphs. 

The way to obtain graphs for 2^3 scattering is by changing the ki line from outgoing into incoming from the 
gluon wave function, as shown in Fig. |6] Naively, one recovers only one graph from the wave functions, which is 
shown in this figure. However, since our calculation is being done in the LFPT, we need to find all possible time- 
orderings. Thus, from the sole graph depicted on the left hand side of Fig. |6]we obtain three graphs of the same 
topology, which are depicted in Fig.]?] This continuation and summation needs to the performed with all the graphs 
in the *I'4 wave function to arrive at all the corresponding graphs needed to calculate M2^3, depicted in Figs. ITlfTOl 
Here we have excluded all the vacuum diagrams since they vanish on the light-front. 

Figs. iTlfTOl group together all the graphs with the same topology. These groups are in turn composed of subgroups 
where the directions of the lines are taken into account. We know that for a graph to have a non-zero value all its 
lines must have a positive longitudinal momentum. Thus, different kinematical conditions which are satisfied by the 
longitudinal momenta of the external lines apply for the different subgroups inside a group. To make this statement 
more clear we can introduce the following classification for the graphs. Let us call G,- the group of graphs which have 
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Figure 6: By changing the outgoing to incoming momentum one can recover the graphs of the amphtudes from the 
wave functions. 




(a) Subgroup Gi_i (b) Subgroup G1.2 

Figure 7: A group of topologically equivalent graphs for 2 — > 3 ampHtude obtained from the continuation of the ki 
momentum, as shown in Fig. |6] 

the same topology, with G2 being the group depicted in Fig. [8] We denote by G, j the subgroup of graphs in group G, 
which are distinguished by different kinematical conditions on the external momenta. In our example G2,; denote the 
subgroups listed in plots a-c in Fig. |8] 
Let 

S2J = sum of graphs in subgroup G2,/ , (51) 
82= sum of graphs in group G2 . (52) 

We can then write 

^2 = S2.l&(Z0 - Z4) @(Z2 - Zl) + S2,2&{Z4 ' Zo) 0(Zl - Z2) + S2,3®(Zo ' Z4) @(Zl - Z2) • (53) 

However, it can be demonstrated that contributions from aU the graphs are actually the same, i.e. ^2 = ■Si.i - ^2,2 - ■S2,3 




(c) Subgroup G2,3 

Figure 8: A group of topologically equivalent graphs for 2^3 amplitude with different time-orderings. 
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(a) Subgroup G3.1 

*1 k-2 




(b) Subgroup G3_2 

Figure 9: A group of topologically equivalent graphs for 2-^3 amplitude with different time-orderings. 





Figure 10: s - channel graph contributions to 2 
respectively. Both consist of only one graph. 



3 scattering. Groups G5 and Gg are shown on the left and right 



(this is discussed in Appendix for a general case). Thus, 

^2 = *2 



(54) 



This sum contains all physical arrangements of the kinematical conditions on the external lines, so we can remove the 
functions. 

We should note that some groups consist of only a single subgroup. These will be composed of only s-channel 
graphs and the internal lines are guaranteed to have a positive longitudinal momentum regardless of the external lines. 
These are graphs shown in Fig. \TQ\ In every group (including the ones with a single subgroup), all of the m graphs 
in one of the subgroups can be obtained from m distinctive graphs belonging to the wave function *1'4. Furthermore, 
if the graphs in one group can be obtained from a certain set of graphs belonging to *1'4, then the graphs in another 
group have to be obtained from a completely different set of graphs. Since the value of all subgroups inside a group 
are the same, in order to calculate M2^„ all we need are the n! graphs that can be obtained by simply changing the 
direction of the ki line in each of the «! graphs in the 1 to « + 1 gluon wave function ^,,+1. We demonstrate this step 
for arbitrary number of gluons below. We shall show that the expression for the value of a graph G is equal to the 
negative value of the expression of the 1 to « + 1 graph H from which G was obtained. In general we know, Ref. lll9ll . 
that a graph wiU be equal to (modulo overall factors) 



ZD ' 



with 



-v. n 



Vk 



z 



n 



(55) 

(56) 
(57) 



internal lines 
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2 to n 1 to n + 1 









- 2j/=i 


A:. 
— ( 


ZO +Z\ 


- Zj,-=2 Zl 


ZO 




Zi 


k- k^ 

ZO Zl 




ZO 


- 2j,=i 


—1 

Zi 



Table 1: Physical conditions for gluon scattering for the case of the 2 — > n and 1 — » n + 1 amplitudes. 

intermediate states 

Thus, we need to compare X, 'V and D for G and H. 

It is important to set up the physical conditions for each of the cases. These conditions are listed in Table[T] Note 
that we have included the energy condition, which we have arrived at by setting the initial and final states on-shell. 

We shall introduce two variables and za- These will be defined differently for the two cases we are working on. 
For the 2 — » « case = -k^ , za = -Zi whereas for 1 — » « + 1 case k^ = k^, Za = Z\- We can now see that these 
definitions allow us to use the same mathematical expression for both cases. To be precise, for both 2 to n and 1 to 
n + I the physical conditions are 

n+l 

i=2 
n+l 

Z0= Za+^ Zi 



^ . £ + y (59) 

Za Za Zi 

Now we can proceed to compare factors Z,- Let's take as an example the graphs in Fig. |6] Looking at these it is 
clear that, using kA and za, we can write, for both cases, 

= ^ = +^2 Zi^Z{)-Za^Za+Z2+Z3, (60) 
kit ^ - !U ~ lU ^ lil !i3^ Zg^ZQ-Z4-ZA^Z2+Z3- (61) 

Once again we see that we can use the same expression for both cases. This is true for any internal line one considers. 
Thus, Z does not change from graphs 2 — > « and 1 ii + I. 

Let's consider the vertex factors 'V. By taking a graph and changing the direction of the ki line only one vertex is 
changed, the one with the ki line (vi). The other vertices remain the same and, because the expressions for the internal 
lines do not change, neither do the expression for these vertices. However, the expression for vi does change: it has 
an extra minus sign, i.e., vi_2^;i = -vi, Therefore the vertex factor contains an extra sign, ^2^11 = -^i^n+i- 

Finally, let us consider the energy denominators £). Let us consider a general 2 — > 2 graph with arbitrary number of 
vertices and internal lines, the graph can also contain loops, it is schematically depicted in Fig. [TT] Let Q denote the 
point where gluon 1 interacts with the rest of a graph. Using it as a reference, we split the graph into a left sub-graph 
and a right sub-graph. Next, let denote by Da and Di, the energy denominators of the intermediate states a and b which 
are present just next to the vertex, we thus assume that there are more vertices to the left and right of the vertex Q, 
such that a and b are the intermediate states. For Fig. Illal we can write 

z),.£„-^.£„ + S_yg.£^_ii_yg, (62) 

ZO Zl Zi Za ^ Zi 
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and 

k- k^ k^ k- ^k^ 

Dt^E,,---^^Eb-- + ^^Eh-y^^, (63) 

zq z\ ZD Za Zi 

where and Eh are the sums of the energies in intermediate states a and b respectively. Here, we've used (|59] | and 
expUcitly excluded the energy of gluon 1 from Ea- Similarly, for Fig. II Ibl we can write 

D„.£,-M.£„-S_yg.£^_i_yS, (64) 

ZO Zl j^Zi ZA j^Zi 

k^ k^ k^ ^ k^ 

Db^Eh-^ + ^^Eb-^ + ^^Eb-)^^ . (65) 

Zq Zl Zo Za Zi 

This time we've excluded the energy of the gluon 1 from Eb- Looking at Da and Db for both cases once again one sees 
that, using fc^ and za, the expressions are the same, though in each of these cases k^ and za are defined differently. 
Hence, the contribution from energy denominators D also does not change. This result is true for all graphs in light- 
front theory with arbitrary topology. We can therefore conclude that the 1 — > n + 1 and 2 ^ n scattering amplitudes 
satisfy the relation 




mn? fen kowwy 



(a) 




mwkn 



(b) 



Figure 1 1 : Schematic representation of the crossing of the momentum of gluon 1 to relate the graphs for 1 — > n + 1 
with graphs 2 — > n. Point Q denotes the vertex at which gluon 1 attaches to the graph. Da and Db denote the energy 
denominators for the intermediate states. 



M2^„({^),zo;^i,zi);{fc2'Z2; ■ • ■ ;^„,z„)) = -Mi^„+i{{k^^,zo}; [k^,ZA,k^,Z2; . . ■ ;^,,z„))k^^-i,,z,,^-z,. (66) 

Furthermore, since the physical conditions for both cases are the same when written in terms of ^ and za, we can use 
them to simplify M2^„ by simplifying M\^„+\ first as much as possible, which we are about to do. 

The usefulness of (|66] | is that it converts a complicated problem into a simpler one (in addition, we get two 
results for the price of one). Generally, M\^n+\ is easier to calculate than M2^„, as we saw previously there are 
smaller number of graphs to compute. Here we will see that it is a relatively simple task to show that M2^„ - 
for the configuration of helicities (-, +,+,■■■ , +) using ( |66] | and the recursion relations derived earlier for the wave 
functions. All one needs to do is to show that Mi^„+i = 0. A substantial step towards this was done in UtIi . and it has 
been summarized in previous section. Note that, unlike the graph shown in Fig. [T] the external lines 1,2, . . . ,n + 1 in 
Mi^„+i will be put on-shell, i.e. there is no energy denominator for the final state. We thus obtain 

n+l 

Mi^„+i(kQ, kuki,..., fe„+i) = -2g ^ ^jzoZlZ2 ■ ■ ■Z(m+i) . ■ ■Z„+ie*"^*v,._i , ^nih,- ■ ■ , /,..., ^„+i) , (67) 

1=2 

17 



where the square root ehminates the extra factors introduced by (|7]i. There is also an overall summation performed 
over all the possible splittings of « gluons of the original wave function. is given by equation ( fT9] l 



^„(\,2,...,n) = (-ir'(/g)''-'A(«> ^ ^ 



-\/Z\Z2 ■ ■ - Zn ^(12...«-l);i^(12...n-2)(«-l;i) ■•■ ^l(2...n) 
1 



. (68) 

V(12...n-1)« V(12...n-2)(«-l n) ■•■ 1'I(2...h-1) 

Groupin g th e terms in (l67l i to the common denominator using the explicit expression ( l68l l one arrives at the following 



formula llTIl (which should be compared with (flST l and (fT9] l ) 
Mi^„iko,ku...,k„)^ (69) 



■ n 2 V^2"=i v*,._^i^(i..ix;+i...«+i)V(i.../x/+i..."+i) 1 



^(12...n)n+l ^(12...n-l)(n«+l) ■ ■ ■ ^1(2.. .n+1) V(12...n)n+1 V(12...n-l)in n+1) ■ ■ ■ Vl(2...«+1) 

(70) 

Using ( l23T l one can show that the sum in the numerator coincides with the difference of the energies 

" I ( " 

The right-hand side of ( fTlT l is proportional to the difference of the incoming light-front energies and the outgoing 
light-front energies, which, according to the physical conditions of the problem, are the same. Hence, the numerator 
in dTOb is equal to zero, which gives Mi^„+i = — > M2^„ ~ Pj^^ - f'j,,, = 0. Note that in arriving to the result for 
the resummed wave function ( |68] | one utilizes the above relation dTTT i for the last intermediate state. In such case the 
sum ( fTTT i is non-zero and it cancels with the denominator giving the non-vanishing resummed relation (|68] |. Therefore 
we see how the recursion relations for the {+;+ + ■ ■ ■ +) wave functions are related to the vanishing of the MHV 
amplitudes with helicities (-, +,+,...,+). 

On the other hand the vanishing of the amplitudes m(+, +,...,+) on the light-front can be argued from the angular 



momentum conservation as follows (see 112511 ) 



The hght-front wave functions satisfy non-trivial conservation of the total angular momentum 

n n—l 



The superscripts z denote the projection along the light-front direction z- The first term corresponds to the contribution 
originating from the spin of the « component Fock state of the wave function in question, whereas the second term is 
the sum over the relative orbital angular momenta. The second term excludes the contribution from the motion of the 
center of mass and is only due to the relative motion. The sum rule ( l72l i has been verified for in QED for a systems 



of electron and spin 1(0) boson [i26[l and in QCD for dressed quark and gluon Il27ll28ll . For a transition of 2 to n the 
relative change of the contribution from intrinsic spin is equal to AS- - n + 2. On the other hand for tree-level g" 
amplitude there are only n triple gluon vertices, each of which can change the angular momentum by Al^ ± 1 . Thus 
the orbital angular momentum can only be changed at most by AL, = n. Thus the amplitude m{+, +,...,+) needs to 
vanish on the light-front. Moreover the vanishing has to occur graph by graph in perturbation theory. 

The same argument for the helicity flip amplitude m{-, +,+,...,+) considered above is not sufficient as the re- 
quirement (|72] | can be in principle satisfied for each graphs. Thus the contributions from individual graphs on the 
Ught-front do not vanish as we saw in the the simple example earlier It is only the sum of graphs which vanishes, and 
the proportionality of the single helicity flip amplitude to the vanishing denominator on the light-front suggests the 



kinematical constraint due to the Lorentz invariance, similarly to the case of the tree amplitudes in gravity 112911 
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3.3. The 2 — > 2 amplitude m{—, +, +) obtained from wave function +, +) by using the recursion relation 

In this section we are going to see how the previous discussion on the relations between wave functions and 
the ampUtudes works in a less trivial case. We shall namely consider the case of (h — > -++) wave function and 
corresponding (+H — >++) amplitude. We shall work with off-shell objects, leaving the comparison with the on-shell 
Parke-Taylor amplitudes to the very end, when the on-shell condition will be imposed. We shall label the wave 
function of the graph in Fig[T]as *P+ + for (1, . . . , «) = (+,..., +) and *I',7+ + for (1, 2, ... , «) -(-,+,...,+). With 
this notation, we can write the amplitude Mi_^„(-i — >- + + ... H-) as 



Mi^„(+ ^ - + + ...+) ^ »P,^:2+(123, 4, ...,«) X / ig-" + ig^ 



loo / \J 

zm \ to - zm) 

+ ^F-f (12, 3, ...,«) X /'"^-•-'" (2g Z2 v(,2)i) 

V Z\2 

+ ^-Yil2,X...,n)x ^l^^^{2gz, vl,,,) 

n-\ 1 

+ Z^"-i"(l'2,...,M-+l,...,«)x J '""-;-'" (2gz,(/.i)n,.,i),) , (73) 
where the terms in parenthesis in the first three lines come from the vertices in Figs. [T3a]-[T3d] respectively. The 
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(c) (d) 

Figure 12: A schematic representation of four terms on the right hand side of recursion (iTJt . The blob in the four 
gluon vertex in graph (a) denotes the combined 4-gluon and Coulomb vertex. Vertical dashed lines indicate that the 
states are taken off'-shell and the energy denominators are taken into account. 

recursion relation (iTJt is schematically shown in Fig. [12] We introduce here notation M to denote that the on-shell 
condition on the last n gluons has not been yet applied. The on-shell amplitudes will be denoted by M. In the case when 
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(a) (b) (c) (d) 

Figure 13: Vertices that are included in recursion relation for the wave functions (l73T l. The gluons on the left hand 
side of the vertices are treated as incoming. 

the energy condition is not applied, then we also can find the gluon wave function *P,7^- from M\^„{-\ — > - + +...+) 
by the following relation 

...,«)= Mi^„(+ ^ - + + ...+) . (74) 

^jZQZ\ ■ ■ ■ z„u„ 

Eqs. ( l73T l and ( f74] i form a recursion relation for the wave functions *1',7^ ^. The first term in Eq. (iTJt contains both 
the 4-gluon vertex and the Coulomb term (which are depicted in Fig. [13]), which have been combined to one effective 
vertex. The recursion relation (iTJt is much more complicated in form then the recursion relation for 'l',^^ ^ written 
previously, compare Eq. |9] The complication stems from the fact that now the recursion involves two distinct wave 
functions for different set of helicities and ^l^'"^ as well as for different values of m due to the presence of 

the 4-gluon vertex and the Coulomb term. This makes searching for the solution to this recursion for general n much 
more involved than in the previous case when the recursion only involved the same type of objects i.e. 

For now let us therefore demonstrate that the general relation derived in the previous section works for the (h — > 
-++) case and that it coincides with the well known results obtained in the literature previously. For this case we have 



Mi^3(+ ^ - + +) ^ n(123) X ^'^ (Z0Z2 - ZIZ3) + %W2, 3) X J^(2gZ2 V(i2)l) 

V 2123 (ZO - Z3)- V Zl2 

+ V(12,3)x ^J^i^gzi v;(i2)) + V(l'23)x ^^[2gz2i v\^) . (75) 
The *P| and ^Pj^ wave functions have already been computed and their expressions are given by Eq. ( fT9] l. The new 



object '^^n be found by multiplying the vertex in Fig. I13cl bv , 



V(1.2)^^-g^^ , \ ■ (76) 

VM21)1 S(21)1^'(21)1 



In addition, we can combine the following terms 

,2 



>y|(123)x ^'^ ZQZ2 + ^r(12.3)x l^{2g Z2v^l2)^) ^%Hl2,2,)x [^{2gz2Voi), (77) 

(Zo - Z3r V Zl2 V Zl2 

^1(123) X ^-^^(-ziZ3) + >I'-+(12, 3) X yj^i^g zi V2(i2)) = >I'r(12, 3) x ^^^(2^ ^i V23) ■ (78) 



Putting in the explicit expressions for "P's the following expression is obtained 



Mi^3(+ ^ - + +) = 2ig^ 



* 2 * 

Z0Z2 VQI _^ ZlZ3 '^23 _|_ ^23 ^32 



Z3Zn V30 zozn v^^ zozi Vqj 



2ig^ 



Z0Z1Z3Z12 V30 



[ziZ2zl VoiV^iV;;, + {ziZ3f vl^V30vli + 223^12^3 v;2V30V;o] , 

^30' 01 



(79) 
(80) 
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where we used Zi2V3(i2) - Z0V30 and Zi2Vq^^2) ~ -^3^30 '■° S^*- '-^^ ^^^^ ^'^^ some manipulations to the 

third term of the second Une: 

^23^1223 Vj2V30vlo = Z23Z12Z3Z2V32V32V30 + Z23Zl2Z3Zlv52Vi5i V30 (81) 

= Z1Z2Z3212 V32 vJjVoi + 22^3^12^32 V32V32 + (ZlZ3)^ vJjVq, V30 + ZoZlZ2Z3 V32 Vqi V30 (82) 
= ZlZ2Z3Zl2vl2vl2^Ql + 22^3ZOV32V3oV32 + 22^3^1 V52V12V32 

+ (ZlZ3)^V32vJiV30 + Z0Z1Z2Z3VI2VQ1V3Q (83) 

where we have used the following relations Z23V30 = Zivoi + Z2V32, ^12^32 = Zov^q + ziv*2 and 223^12 - ZiZj, + Z0Z2 in 
order to rewrite the expressions. Thus, 

Mi^ii+ ^ - + +) = ^-r[ZlZ2Z3Zl2vl2^l2^01 + (ziZ2zl V01V01V30 + 22^3^0^32 V30V32) 

ZQZiZiZn ^3oV3oVqj ^ ' 

+ (222321^52^12^32 + 20212223 V32V()iV3o)] (84) 

Finally, we can group some terms which are proportional to denominator D3 

2l222o mVoi^lo + 22232oV32V'3oV32 = -^2o22223V'30^3, (85) 

2 1 

222321 V32V12V32 + 20212223V32V01V30 = --212223V3203- (86) 



Therefore, 



Mi^3(+ ^ - + +) - ^ 



202l232l2 ^'30V3oVqj 



212223212V32V52V01 - ^Z2D3 (2123 v52 + 20223 V30) 



(87) 



The on-shell energy condition D3 = can now be applied to the above expression as required for the on-shell 
amplitudes. This leaves only the first term in Eq. (l87t which is finite and it matches the Parke-Taylor amplitude. This 
can be made more explicit by utiUzing the relation 2oVgiV3o = -22VJ2V32 and 2o2iVqjVoi = 2223V23V23 which leads us to 
the expression 



2ig\zoZi) 



Mi^3(+ ^ - + +) = -2.^1:^ 01 (88) 

202 12223 V01V12V23V30 

which is the exact result for this case. Thus by utilizing the recursion relation for the wave functions we are able to 
reconstruct the the on-shell Parke-Taylor amplitude. The full proof would of course involve solving exactly the recur- 
sion relation, which we leave for a future work. An interesting structure though emerges from this example, namely 
the fact that the ofF-shell amplitude M can be written as a sum of the on-shell amplitude plus the term proportional 
to the energy denominator which vanishes if one takes the on-shell condition D3 — > 0. Thus we can expect that the 
general structure for arbitrary number of final particles in the LFPT has the following form 

M(2 «) = M(2 ^ «) + 0(D„) , (89) 

where the term which is related to the off-shellness of the amplitude vanishes upon taking the physical condition 
Pt = Pi.. , i.e. it is at least linear in the energy denominator D„. 



4. Conclusions and outlook 

In this paper we have investigated the gluon wave functions, fragmentation functions and scattering amplitudes 
within the framework of the light front perturbation theory. The recursion relations for each of these objects have 
been constructed on the light-front. For the special case, when the helicities of the outgoing gluons are the same, it is 
possible to solve the recursion relations for the gluon wave functions and fragmentation functions and obtain the final 
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results in a very compact form. Furthermore, we have shown that in general the scattering amplitudes (with ofF-shell 
final states) can be computed from the gluon wave functions. One of the advantages of using the LFPT is the fact that 
the variables used to express the helicity amplitudes naturally arise in this framework. 

In particular, we have shown that the compact recursion relations for the gluon wave functions derived before for 
the special case of the helicity configuration are related to the vanishing property of the on-shell helicity amplitudes 
for selected configurations of the helicities. We have generalized the recursion relations for the wave functions and 
the fragmentation functions to include general configurations of helicities. In the latter case, the recursion relations 
are light-front analogs of the previously used Berends-Giele recursion relations. We have also verified the vanishing 
of the amplitudes on the light-front for special cases of the helicities, which originates from the angular momentum 
conservation and energy conservation. Using the general relations between and gluon wave functions and scattering 
amplitudes we have been also able to reproduce some of the results available in the literature. The light-front methods 
presented here can be used to compute the scattering amplitudes with arbitrary number of external legs and for di fferent 
helicity configurations. 
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Appendix 




(a) (b) 

Figure 14: Generic 2 to n tree-level graphs. In (a) gluon c departs from the lower subgraph at Q and arrives at the 
upper subgraph at P. For (b) the opposite is true. 

In this Appendix we shall investigate a general structure of the tree-level amplitudes in light-front theory. In par- 
ticular we shall demonstrate that the typical tree-level graphs can be written in terms of the factorizable contributions 
originating from initial and final state emissions and non-factorizable contribution which typically comes from the 
energy denominator of the intermediate state exchanged between two subsets of graphs. On top of that we shall show 
for a given group of topologically identical graphs, all the subgroups which differ by the orientation of the exchanged 
particle within such group have the same value. This is true with with the exception of those subgraphs for which 
the physical conditions constrain the graphs to have a value of zero (i.e. vacuum graphs). We will show this by 
demonstrating that graphs within the given subgroup sum up to the expression which can be analytically continued to 
obtain the graphs from the different subgroup. Generic tree-level graphs are shown in Fig. I14al and Fig. ll4bl In these 
figures, the top and bottom box both have a set of vertices and thus form subgraphs with a definite topology. The 
graph represents the sum of all the possible time-orderings between the upper and lower parts of the diagrams. 

Points P and Q denote the vertices at which gluon labeled by c attaches to both subgraphs. We need also to perform 
the sum over the relative orderings of these vertices with respect to the orderings in the lower and upper subgraphs. 
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For instance, for situation depicted in Fig ll4al we need to perform the sum over all the orderings of the splittings in 
the upper box which occur before the vertex P with respect to the position of vertex Q. One can call them initial state 
splittings of the subgraph U, as they occur before the interaction P. Similarly, we have to sum over all the orderings 
of the splittings of lower subgraph L with respect to P which occur after interaction depicted by the vertex Q. These 
splittings one can refer to as the final state splittings of the subgraph L. A similar situation is valid for a graph shown 
inFig.[T4bl 




Figure 15: 2 to n tree-level graph with only final state splittings. 



Let us consider a graph with only final state splittings, that is with the splittings which occur in the upper part 
occur only after vertex P and splittings in the lower subgraph occur only after vertex Q. Such graph is depicted in 
Fig. [15] Let us label by A, the energy denominators of the upper subgraph U in the absence of lower subgraph and 
gluon c, and by Bj the energy denominators in the lower subgraph in the absence of the upper part as well as the 
gluon c. In particular A,(fiy) is the energy denominator of the intermediate state which is on the position when we 
count from the right hand side of the corresponding subgraph. We denote by C the energy of the intermediate gluon 
c. Consider a class of subset of graphs in which all the splittings in L and U occur after vertex P. It was shown in 
1 1711 that in such case the fragmentation trees factorize. Using this property for such configuration we can write the 
contribution from the energy denominators in the form 



1 



1 



C + BAiA2...A„,Bi 



■ B„ 



(90) 



where by B we denoted the difference of the energy denominator of the gluon 0,i.e. B - Bq — Bq where Bq is the 
energy of the incoming gluon into the lower subgraph L and Bq is the sum of the energies of the final states in the 
same subgraph L. 

The next step is to include the remaining graphs of the same topology, but this time with the splittings in the 
lower subgraphs occurring earlier than the point P but later than Q. Summing this subset of diagrams together with 
the diagrams which result in ( |9Ut one arrives at the factorized expression. Again, the upper tree factorizes from the 
lower part and the sum of all the time-orderings for the graphs in Fig. [TSlis 



1 



1 



-C +BA^A2...A„,Bi 



■ B„ 



The first factor is just the energy denominator for the simple 2 — » 2 graph, without any splittings in the final 
state, with the only assignment that the final energy of the lower leg in this graph is given by Bo- This denominator in 
turn multiplies the factorized contributions to the fragmentation trees originating from subgraphs L and U. The same 
procedure can be repeated for the general graph depicted in Fig.[T6]with only the initial state radiation. Finally one 
needs to consider the most general case in which the splittings can occur both before and after P and Q on both the 
upper and lower subgraphs simultaneously. In this case the contribution from the energy denominators is given by 



1 



1 



1 



A1A2 . . .AkBiBi ...Bi-C + B A1A2 . . .A,„Bi 



(91) 
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Figure 16: 2 to n tree-level graph with only initial state splittings. 



A similar result can be obtained for the graph Fig. II 4a I f or which the contribution from the energy denominators reads 

1 1 1 

A1A2 . . .AkBiB2 ...B,-C+A A1A2 . . .A,„Bi ...B„ ' 

where now A is the difference between the energy of the incoming gluon 1 in the upper subgraph U and the sum of 
the final energies of the states in the same subgraph U. Using global light-front energy conservation A = -B and 
performing substitution C - -C we see that the contribution from the energy denominators for both graphs is the 
same up to a global sign. One needs to take into account also the additional sign change between graphs [T4bl and [T4al 
coming from the contributions from the longitudinal fractions. Thus the general expression for the graph (either [l"4a1 
or ll4bl i has a compact (partially) factorizable form with the energy denominators given by expression (|911 . 

References 

[1] S. J. Parke and T. R. Taylor, Phys. Rev. Lett. 56, 2459 (1986). 

[2] F. A. Berends and W. T. Giele, Nucl. Phys. B 306, 759 (1988). 

[3] M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau and A. D. Polosa, JHEP 0307, 001 (2003) |hep-ph/0206"293l . 

[4] T. Gleisberg and S. Hoeche, JHEP 0812, 039 (2008) | arXiv:0808.3674 [hep-ph]]. 

[5] R. Kleiss and G. van den Oord, Comput. Phys. Commun. 182, 435 (201 1) | ar Xiv: 1006.5614] [hep-ph]]. 

[6] A. Cafarella, C. G. Papadopoulos and M. Worek, Comput. Phys. Commun. 180, 1941 (2009) larXiv: 07l0.2427] rhep-ph11. 

[7] L. J. Dixon, J. M. Henn, J. Plefka and T. Schuster, JHEP 1101, 035 (201 1) |arXiv:1010.3991 [hep- ph]]. 

[8] J. Alwall, M. Herquet, R Maltoni, O. Mattelaer and T. Stelzer, JHEP 1106, 128 (201 1) I arXiv: 1 1 06^0522] [hep-ph] 1 . 

[9] R A. M. Dirac, Rev. Mod. Phys. 21, 392 (1949). 

[10] L. Sus.skind, Phys. Rev. 165, 1535 (1968). 

[11] J. B. Kogut and L. Susskind. Phys. Rept. 8, 75 (1973). 

[12] S. Weinberg, Phys. Rev. 150, 1313 (1966). 

[13] S. J. Brodsky. C. D. Roberts, R. Shrock and P C. Tandy, Phys. Rev. C 82, 022201 (20 10) larXiv:1005.4610| [nucl-th]]. 

[14] S. J. Brodsky and G. F de Teramond, Phys. Lett. B 582, 21 1 (2004) |hep-th/03 102271 - 

[15] A. H. Mueller, Nucl. Phys. B 415, 373 (1994). 

[16] R Antonuccio, S. J. Brodsky and S. Dalley, Phys. Lett. B 412, 104 (1997) |hep-ph/97054T3]. 

[17] L. Motyka and A. M. Stasto, Phys. Rev D 79, 085016 (2009) [arXiv:hep-ph/090 1.4949]. 

[18] J. B. Kogut and D. E. Soper, Phys. Rev. D 1, 2901 (1970). 

[19] J. D. Bjorken, J. B. Kogut and D. E. Soper, Phys. Rev. D 3, 1382 (1971). 

[20] G. P Lepage and S. J. Brodsky, Phys. Rev. D 22, 2157 (1980). 

[21] S. J. Brodsky, H. C. PauH and S. S. Pinsky, Phys. Rept. 301, 299 (1998) rarXiv:hep-ph/97054771. 

[22] M. L. Mangano and S. J. Pai'ke, Phys. Rept. 200 (1991) 301 | arXiv:hep-th/0509223 |. 

[23] D. A. Kosower, Nucl. Phys. B 335, 23 (1990). 

[24] M. T. Grisara and H. N. Pendleton, Nucl. Phys. B 124, 81 (1977). 

[25] S. J. Brodsky and G. F de Teramond, PoS LC 2010, 070 (2010) | arXiv: 1 1 0.4962 '[hep-th]] . 

[26] S. J. Brodsky, D. S. Hwang, B. -Q. Ma and I. Schmidt, Nucl. Phys. B 593, 311 (2001) | hep-th/0003082| . 

[27] A. Harindranath and R. Kundu, Phys. Rev. D 59, 116013 (1999) |hep-ph/9802406|. 



24 



[28] B. -Q. Ma and I. Schmidt, Phys. Rev. D 58, 096008 (1998) jhep-ph/9808202l . 
[29] M. T. Grisara, P. van Nieuwenhuizen and C. C. Wu, Phys. RevTD 12, 397 (1975). 



25 



